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We analyze the dependence of heavy quark free energies on the baryon chemical potential fit 
in 2-flavour QCD using improved (p4) staggered fermions with a bare quark mass of m/T — 0.4. 
By performing a 6 th order Taylor expansion in the chemical potential which circumvents the sign 
problem. The Taylor expansion coefficients of colour singlet and colour averaged free energies 
are calculated and from this the expansion coefficients for the corresponding screening masses are 
determined. We find that for small fib the free energies of a static quark anti-quark pair decrease 
in a medium with a net excess of quarks and that screening is well described by a screening mass 
which increases with increasing fib. The /^-dependent corrections to the screening masses are well 
described by perturbation theory for T>2T C . In particular, we find for all temperatures above T c 
that the expansion coefficients for singlet and colour averaged screening masses differ by a factor 2. 

PACS numbers: 11.15.Ha, ll.10.Wx, 12.38Gc, 12.38.Mh 

I. INTRODUCTION 

Numerical studies of QCD provided quite detailed information about the properties of matter at high temperature 
and vanishing net baryon density [1]. In particular, the screening of static quark anti-quark sources at large distances 
and their renormalization has been analyzed in quite some detail [2-4]. Compared to this our knowledge on the 
dependence of the equation of state and screening at non-zero baryon number density, or equivalently, non-zero 
baryon chemical potential (fib) is rather limited. The /^-dependence of the QCD partition function [5] and the HTL- 
resummation of the pressure [6] have been evaluated only recently. Although in leading order of high temperature 
perturbation theory the dependence of the Debye screening mass on /if, is well-known [7] neither the temperature 
range for the validity of this perturbative result nor the generic features of screening of heavy quark free energies at 
non-zero fib have been analyzed so far with non-perturbative methods in the vicinity of the QCD transition . 

Recently studies of the equation of state have successfully been extended to non- vanishing baryon chemical potential 
using Taylor expansions [9] around fib = as well as reweighting techniques [10] and imaginary chemical potentials 
[11]. We will use here the former approach to analyze the screening of static quark anti-quark sources at non-zero fib, 
i.e. in a medium with a non- vanishing net quark density. We evaluate the Taylor expansion coefficients for correlation 
functions of heavy quark anti-quark pairs and deduce from this expansion coefficients for the screening mass. 

For this work we analyzed gauge field configurations using a p4-improved staggered fermion action with Nf = 2 
degenerate quark flavours. We used the same data sample that recently has been generated by the Bielefeld-Swansea 
collaboration for the analysis of the equation of state [9]. This sample consists of 1000 up to 4000 gauge field 
configurations available for several bare gauge couplings below and above the transition temperature. The lattice 
size is 16 3 x 4 and the bare quark mass, rh/T — 0.4, corresponds to a pion mass of about 770 MeV. In addition we 
generated 1000 configurations at T = 3T C and 1600 at T — 4T C to check for the approach to the high temperature 
perturbative regime. In addition to gauge invariant colour averaged free energies we also have analyzed colour singlet 
free energies. To do so all gauge configurations have been transformed to Coulomb gauge before evaluating the 
Polyakov loop correlation functions. 

This paper is organized as follows. In the next section we discuss the general setup for calculating Taylor expansions 
for heavy quark free energies. Our numerical results for singlet and colour averaged free energies are presented in 
section III. In section IV we discuss the determination of //^-dependent corrections to screening masses from the free 
energies. Our conclusions are given in section V. In an appendix we give detailed expressions for the Taylor expansion 
coefficients of purely gluonic observables. 



A first attempt to calculate heavy quark free energies at non-zero quark chemical potential has been discussed in [8]. 
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FIG. 1: The O" 1 order coefScients fgQ and /q 1 ^ for the singlet and colour averaged free energy in the vicinity of T c . Jqq 
is matched to the T = heavy quark potential at small distances (a). 



II. TAYLOR EXPANSION OF HEAVY QUARK FREE ENERGIES 



A heavy (static) quark Q at site x is represented by the Polyakov loop, 



L(x) = Y\ U 4 (x 7 x 4 ) 



(1) 



X4—1 



which is an SU(3) matrix. A heavy anti-quark Q is described by the corresponding hermitian conjugate matrix. 
The free energy of a QQ-pair separated by a distance r is then calculated from the expectation value of the cor- 
relation function of L(0) and L^(r) where r points to a site with distance r to 0. The dependence on the baryon 
chemical potential fib or quark chemical potential /i = /if,/3 is established solely through the fermion determinant, 
det M ({Up(x)}, fi, rh), with in denoting the bare quark mass. In order to avoid the sign problem, which arises from 
ite(detM) not being positive definite for fi ^ 0, Taylor expansions in the quark chemical potential are used. This 
allows us to perform our simulations at zero chemical potential thereby restricting us to small chemical potentials. 

A purely gluonic observable O like the Polyakov loop L(x) or a corresponding correlation function does not explicitly 
depend on the quark chemical potential; it is calculated in terms of link variables U p {x) of the gauge field configuration 
which do not explicitly depend on fi. Any /i-dependence of the expectation value (O) thus arises from the fx- 
dcpcndcncc of the Boltzmann weights in the QCD partition function, i.e. the //-dependence of the fermion determinant. 
Therefore we can apply the same method that was used for the power series expansion of the equation of state; 
expanding the fermion determinant in powers of /i leads to a power series of our purely gluonic observable, 



(O) =(0>-(l + 0l/ i + o 2At 2 + .-.) 



(2) 



where (O) denotes the expectation value of O evaluated for vanishing chemical potential. We consider observables 
like the colour averaged and singlet QQ-correlation functions, 



O av (r) 



01 M = TrW y E TrL to Li to> 



Af A, 



(3) 



where the sum refers to all sites x, y with \\x — y\\ = r and N is the number of these x, y-pairs. As O av ^ and the 
corresponding expectation values are strictly real for every single gauge field configuration the odd orders in the 
expansion vanish as was argued in [12]. For observables like the Polyakov loop itself or static quark-quark correlations 
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FIG. 2: The 2 nd order coefficients of the singlet (a,c) and colour averaged (b,d) free energies for some selected temperatures 
below (a,b) and above (c,d) T c . 



like TrL(0)TrL(r) we also have to take into account the odd orders which are in general non- vanishing. In the 
appendix we give explicit formulas for calculating the expansion coefficients of an arbitrary gluonic observable up to 
sixth order in fJ./T. We have used this to evaluate the first three, non-vanishing expansion coefficients of the purely 
real observables considered here. 

We extract the colour averaged free energy of a static quark anti-quark pair from the Polyakov loop correlation 
function 

F™ Q (r,T,ri = -T\n(0™(r))^ (4) 

and the colour singlet free energies from 

F^(r,T )M ) = -Tln(0 1 (r)) /i . (5) 

We renormalize the Polyakov loop as described in [13] such that at short distances and vanishing chemical potential 
the singlet free energy F^q(v, T, 0) matches the zero temperature heavy quark potential. This also fixes the renormal- 
ization of the Polyakov loop and all its correlation functions. In particular, this renormalizes also the colour averaged 
free energies. As all our calculations have been performed on lattices with temporal extent N T — 4 the smallest 
available distance at which this matching could be performed is r$ = 1/4T. 

In order to determine the expansion coefficients of the colour averaged (av) and singlet (1) free energies, 

*3flkT,M) = fQQ^,T) + r QQ , 2 (r,T) ^ + f^(r,T) (£)V^ >6 (r,T) (^f + O [[^ , (6) 
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FIG. 4: The 6 th order coefficients of the singlet (a) and colour averaged (b) free energies for temperatures above T c . 



with x = av and 1, we apply (A18) to the corresponding Polyakov loop correlation functions. We again note that 
these are strictly real on every gauge field configuration and thus have an expansion in even powers of n/T. Explicit 
formulas used for the calculation of the expansion coefficients (r, T) and (r, T) are given in the appendix. 



III. NUMERICAL RESULTS ON QQ FREE ENERGIES 

In Fig. 1- 4 we show the leading and higher order expansion coefficients up to sixth order in fi/T expressed in units 
of the string tension 2 . We do not include data for all temperature values analyzed by us because for T 3> T c they 
have very small absolute values and for T < T c they suffer from large statistical errors and are still consistent with 



2 We use /3 C = 3.649 [9, 12] as the bare pseudo-critical coupling. Results on the string tension [14], ay/a, are then used to set the 
temperature scale in units of T. 
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FIG. 5: The coefficients for the singlet and colour averaged free energies at infinite distance rT versus temperature. They have 
been obtained from a weighted average of fgg n (r,T) and fgg n (r,T) at the five largest distances. 



zero. 

The leading order results, J^q Q (r, T) and Sqq {t, T) are consistent with previous analyses of static quark anti-quark 
free energies performed in 2-flavour QCD at fi = on the same data set [3] . For the second order expansion coefficients 
we display separately results below (Fig. 2(a), (b)) and above (Fig. 2(c), (d)) the /i = transition temperature, T c . 
As can be seen the second order expansion coefficients are always negative and increase in magnitude in the vicinity 
ofT c . 

The corresponding results for the 4 th and 6 th order expansion coefficients are shown in Fig. 3 and Fig. 4, respectively. 
Here we only show results above T c ; below T c the expansion coefficients are consistent with being zero within errors 
even at rather short distances and errors grow large for rT > I. We note that all expansion coefficients shown in 
Figs. 2 to 4 vanish at small distances. This shows that a quark anti-quark pair is not affected by the surrounding 
medium if its size becomes small. This observation also justifies our procedure to renormalize the Polyakov loop 
by matching the [i — singlet free energy to the T = heavy quark potential. The renormalization constant is 
independent of fi. 

Also close to T c , where the /i-dependcncc of the free energies is strongest, the absolute values of the fourth and sixth 
order expansion coefficients are of the same order as or smaller than the second order expansion coefficient. Therefore 
the 4*' 1 and 6 th order contributions rapidly become negligible for [ijT < 1. 

Although the errors are large for the higher order expansion coefficients they show that at high temperature the 
2 nd and 4 t h order expansion coefficients are opposite in sign, fQQ 2 i r ^) < anc ^ $QQ J* 1 "^) > 0- This is consistent 
with the expectation that at high temperature the asymptotic large distance value of the heavy quark free energy 
is proportional to the value of the Debye mass [15]. In this limit one obtains alternating signs of the expansion 
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FIG. 6: The singlet free energies Fqq as function of distance for finite chemical potential and for various temperatures 
coefficients of the heavy quark free energies when one expands the perturbative Debye mass [7], 



m D (T,n) _ m D . (T) 3N f /^\ 2 

g{T)T ~ g{T)T V + {2N C + N f )^ \t) ' 
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with mj) :0 (T) = g(T)Ty ^ + denoting the Dcbye mass for vanishing baryon chemical potential. Although 
the statistical significance of our results for /qq 6 {r, T) rapidly drops with increasing temperature this pattern of 
alternating signs seems to be valid also at sixth order at least for temperatures T>1.05T C . 

Except for temperatures close to the transition temperature the asymptotic behaviour of the free energies is reached 
at distances rT> 1.5. We determined their large distance value by taking the weighted average of the values at the 
five largest distances. The results are shown in Fig. 5. We note that l/gg^l 00 ! ^)l have a pronounced peak at T c . 
This also holds for |/g^ 2 ( r ' evaluated at any fixed distance r. In fact, /q^^ ( r > ^) i s proportional to the second 
derivative of a partition function including a pair of static sources, QQ. It thus shows the characteristic properties of 
a susceptibility in the vicinity of a (phase) transition point. 

Fig. 5 also shows that at large distances, within the statistical errors of our analysis, the expansion coefficients for 
the colour averaged and singlet free energies approach identical values, 

/a V Q,n(°°' T ) = /«kn(^ T )> (8) 

where 

t ) = }™ fQQJ r > T ) ' with x = av ' 1 • ( 9 ) 

This has been noted before at [i = and suggests that at large distances, e.g. for rT>1.5, the quark anti-quark 
sources are screened independently from each other; their relative colour orientation thus becomes irrelevant. 

Including all terms up to sixth order we calculated the singlet and colour averaged free energies in the range from 
[ijT = 0.0 up to 0.8. Results for the colour singlet free energies evaluated at a few values of temperature are shown 
in Fig. 6. Similar results hold for the colour averaged free energies. The free energies decrease relative to their values 
at [ijT = for all temperatures above and below T c . At small distances the curves always agree within errors. With 
increasing distance a gap opens up which reflects the decrease in free energy at non zero /z. As indicated by the 
asymptotic values 2 (oo,T), which give the dominant /i-dependent contribution at large distances, the medium 
effects are largest close to the transition temperature and become smaller with increasing temperature. 



IV. SCREENING MASSES 



For temperatures above T c and large distances r the heavy quark free energies are expected to be screened, 



AF^(r,7» = F%\oo,T,VL)-F$(r,T tli ), 



(10) 



with n = 1, 2 for the singlet and colour averaged free energies respectively. In the infinite distance limit we thus can 
extract the screening masses, 



m 



^(T,^) = - r lim o iln(AF- 1 (r,T, M )) 



(11) 



We use this as our starting point to derive a Taylor expansion for the screening masses. Expanding the logarithm in 
eq. 11 in powers of fi/T it is obvious that also the screening masses are even functions in [i/T, 



»'" v; ■:/•„:. = nC Hi • m^{T)(^) +,„■ 



(12) 



To analyze the approach of the various expansion coefficients to the large distance limits we introduce effective masses, 

TO eff,n( r > T )> with X = &V,1, 



< ff , 2 (r,T) = 
m x eS , 6 (r,T) = 



^foo.^T) 1 f A f X QQ,2(r,T) 



QQA 



A f X QQ^ T ) 



+ 



3 V A W r ' T ) 



(13a) 
(13b) 

(13c) 



8 




2 
1.5 
1 

0.5 


-0.5 
-1 

-1.5 
-2 



T/T =1.11 



-f — i i{ \ I 



rf 

0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2 

(b) 



FIG. 7: Example for the distance dependent effective masses converging at large distance to the expansion coefficients for the 
screening mass which are represented by the horizontal lines. 



In the limit of large distances these relations define the expansion coefficients of the colour averaged and singlet 
screening masses, 

mr 1 (T)= limm^r.T). (14) 

r — >oc ' 

As will become obvious in the following the effective masses defined above show only little r-dependence. They are 
thus suitable for a determination of the /x-dependent corrections to the screening masses. This is not the case for 
the leading order, /i- independent, contribution. In order to determine mJ(T) we use an ansatz for the large distance 
behaviour of the singlet free energy motivated by leading order high temperature perturbation theory, 

W',T) - W°o,T) " |^e-™iP> (15) 

We fit our data to this equation using ao(T) and mJ(T) as fit parameters where /gg (oo, T) is determined as described 
in the previous section. We choose the same fitting procedure as in [3] namely averaging results received from five fit 
windows with left borders between rT = 0.8 and rT = 1.0 and right border at rT = 1.73. While the above ansatz is 
known to describe rather well the large distance behaviour of the color singlet free energy, it also is known that the 
sub-leading power-like corrections are much more difficult to control in the case of the colour averaged free energy. 
For this reason we will analyze here only the leading order contribution to the singlet screening mass. 

Results for effective masses in the singlet channel are shown in Fig. 7 as function of rT for one value of the 
temperature. As can be seen the asymptotic value is indeed reached quickly before the errors grow large at distances 
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FIG. 8: Expansion coefficients of screening masses in the colour singlet (m„) and colour averaged (m^ v ) channels. Except for 
7i = m^T) is determined from the r — > oo limit of (13). The lines are the first order perturbative predictions according to 
eq. 20. The dotted lines in (a) show the la-range of the x 2_r h with parameter A. 

rT>l. The expansion coefficients m^ /,1 (T), m^' (T) and mg V,1 (T) are thus well determined from the plateau values 
of these ratios. Similar results hold in the colour averaged channel. We found the left border of the plateau to lie 
between rT = 0.48 close to T c and rT = 0.23 for T > 1.15T C . Results for the various expansion coefficients are 
shown in Fig. 8. This figure shows that at high temperatures the /^-dependent corrections to the screening mass of 
the colour averaged free energies m av (T, /z) are twice as large as those of the (Debye) screening mass in the singlet 
channel, to 1 (T, /i) . This is expected from perturbation theory, which suggests that the leading order contribution to 
the colour singlet free energy is given by one gluon exchange while the colour averaged free energy is dominated by 
two gluon exchange. Using resummed gluon propagators then leads to screening masses that differ by a factor of 2, 

m x n (T) = lm^(T), n = 2,4,6, (16) 

Our results suggest that this relation holds already close to T c (Fig. 8). We thus have no evidence for large contributions 
from the magnetic sector, which is expected to dominate the screening in the colour averaged channel at asymptotically 
large temperatures [16] and which would violate the simple relation given in eq. 16. 

In order to compare the expansion coefficients with perturbation theory we need to specify the running coupling g(T). 
Following [3] we use the next-to-leading order perturbative result for the running of the coupling with temperature 
but allow for a free overall scale factor. We thus fit our data on the T-dependence of the leading order (/x = 0) 
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2 nd order expansion coefficients in \x (b). The line in (a) shows the leading order perturbative prediction, which is 3/8-7T 2 . 



screening mass by the ansatz, 



m 1 Q (T) = A-^=g(T)T, (17) 



with the 2nd order perturbative running coupling, 

g(T) 2 



V3" 



29 /JTX 115 1 , (V T 



247T 2 VA^; 232^ V 



(18) 



where we use T c /Ajjg — 0.77(21) and the scale ju = 2tt as in [3]. Fitting our data to eq. 17 with fit parameter A, 
yields 

A = 1.397(18) , (19a) 

which is almost identical to the result in [3] where the data for T = 3T C were still missing. Our fit result is included 
in Fig. 8. We also compare the temperature dependence of m\{T), m\(T) and m\(T) with corresponding expansion 
coefficients of the perturbative Debye mass which result from an expansion of (7) using (17) as the 0th order. These 
expansion coefficients are alternating in sign, 



m D ,e(T) = ££e-MT). 



(20a) 
(20b) 
(20c) 



At least for the second order coefficient m\{T) we find that this yields a satisfactory description of the numerical 
results for T>2T C . Eq. 20 shows that subsequent terms differ by about an order of magnitude, which explains why 
our signal for a non-zero contribution m^(T) is rather poor for n > 2. 

From (20) we find mD,2(T)/mn,o(T) = 3/87T 2 which is independent of A and g(T) and is compared with our 
numerical results in Fig. 9(a). We note that the perturbative value for this ratio is already reached for T/T c >2. In 
Fig. 9(b) we show the /i-dependence of the singlet screening mass for a small values of n/T. Here we included only 
contributions from the 0th and 2nd order expansion in the calculation of m 1 (/i, T)/T. 
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V. 



CONCLUSIONS 



We have analyzed the response of colour singlet and colour averaged heavy quark free energies to a non-vanishing 
baryon chemical potential and have calculated the resulting dependence of screening masses on the chemical potential. 
Using a Taylor expansion in fi/T we get stable results for the leading, non- vanishing correction, m\{T), which is 
C((/i/T) 2 ). We find that this correction in absolute units as well as its ratio with the leading order screening mass, 
mJ(T), is large in the vicinity of the transition temperature. The ratio m\{T) j 'mJ(T) is in agreement with perturbation 
theory for T>2T C indicating that the expansion coefficients (T) receive the same multiplicative rescaling as the 
leading order screening mass. 

A calculation of the ^-dependent corrections to the screening mass in the colour averaged channel shows that these 
corrections are twice as large as those in the color singlet channel for all temperatures T > T c . This agreement 
with leading order perturbation theory indeed is quite remarkable as it suggests that the leading contribution to the 
^-dependent corrections of the colour averaged screening mass is due to two-gluon exchange. 

The higher order expansion coefficients of the screening mass vanish within statistical errors at temperatures larger 
than 1.2T C . The analysis of the asymptotic behaviour of the free energies themselves, however, suggests that these 
corrections are non-zero but small at high temperature and have alternating signs. This is consistent with the leading 
order perturbative result for the Debye mass subsequent expansion coefficients of which drop by more than an order 
of magnitude and alternate in sign as they arise from an expansion of a square root. 

Our results thus suggest that at least for small values of the chemical potential and fixed temperature the screening 
length in a baryon rich quark gluon plasma decreases with increasing value of the chemical potential. This is consistent 
with the expectation that the transition to the high temperature phase shifts to lower temperatures at non-zero baryon 
chemical potential. 
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APPENDIX A: CALCULATION OF EXPANSION COEFFICIENTS 



The ^-dependent expectation value of a complex quantity O is 




(Al) 



where is the partition function for finite /i and where A = (det M(p)) sl is the determinant of the fermion 
matrix. In the following we denote expectation values for vanishing /x as (•••) = (••• ) . We define 



(Tin A 



Nf d n In det M 
~4 



fj,=a 



(A2) 
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The L n can be written as traces over the inverse of the fermion matrix and its derivatives 
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\, l 9M »r 1 5M 



dn 3 

on 



20M-^M~^M 
dn 

M 

ctyi 2 



9yU 



j9 3 M 



ju=0 



\d 2 M ^ i<9 4 M !<9 3 Af ,d 3 M 

— 15M M 10M M 

dn dn 5 dn 2 dn 4 dn 3 dn 3 

,9M,, l( 9 4 M ,9M„ ,9 2 M„ , d 3 M id 2 Af ,9M 
+30M — — M — — M — — j- + 60M — — M — — — M -— — - + 60M -—-M — — 

9/U a^i a/i 4 a^i ctyz' 3 a^ 2 on 

d 3 M 



M" 1 - 



— 180M — — M — — 

on on 

+360M" 1 ^— M _1 — — 



Af - — — -M~ — — - 



M" 1 



9yU 2 

M -idM 

dn dn 

— — M -1 ^— M _1 — - 
a^i o^i O/U / 



on 



id 2 M^ ,9M„ ,9 2 M 
_ 90M" 1 — M- 1 —^ M- 1 — M- 1 ^- 

(J 11^ rin ri 1 1 ' 



M -id 2 M 

dn 2 



dn dn 2 dn 

,9M„, ,9M 
120M — — M~ — — 



i^M,^ ,0M W ,9M„ ,a 3 M 
— -5- - 120M" 1 — M' 1 — M- 1 — M- 1 — T - 

^ M _1 - 



1^ 2 

30M~ — -^-M _ — 



9^ 2 

M -5- 
dn 



From 



M\n) = j 5 M{-nh 5 



it follows that L n is real for even and imaginary for odd n. Using A = e lnA we find 



where 



A(n) - A(O)(l + £ 1 /i + £ 2M 2 + --- + A;Ai 6 + 0(M 7 )) , 



D\ = Li , 
D2 = \ {L\ + L 2 ) , 

D 3 - i (L 3 + 3LiL 2 + L 3 
6 v 

l(L 4 + 6i 2 i 2 + 3L 2 



^ 4 = -1 



£>5 = 
^6 = 



, 2 + 4Lii 3 + L 4 ) , 
^ (if + 10L 3 L 2 + YoL\L\ + 10L 2 L 3 + 10L 2j L 3 + 5iiL 4 + £ B ) , 
^ (if + 15L 4 L 2 + 45L 2 L 2 + 15^ + 20L 3 L 3 + &QL X L 2 L 3 
+ 10L 2 + 15L 2 L 4 + 15L 2 L 4 + QLxL 5 + L 6 ) . 



(A7) 



(A8) 

(A9) 

(A10) 

(Alia) 
(Allb) 

(Allc) 

(Alld) 

(Alle) 

(Allf) 
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We immediately see that D n is real for even and imaginary for odd n. Because 

Z (I = (l + % + ...% 6 )+% 7 ) (A12) 

is real one has (D n ) = for odd n. We consider the case where the observable O is independent of \i. The expectation 
value (Al) then becomes 

{0) » 1 + (D 2 ) fi 2 -\ + (D 6 ) + ° {P ' ■ ( 3) 

Expanding in powers of \i we get 

(O)^ = (O) (1 + OiH + (-V 2 + 2 ) f + (-^Oi + 3 ) M 3 + - Z>4 - W2 + 4 ) ^ + (X>§Oi - 

-V 2 3 + 5 ) ^ + (-X>! + 2P 2 I?4 - X> 6 + ^2 - Z> 4 2 - V 2 4 + C 6 ) + 0(fi 7 ) (A14) 

where we use the notation 

Oi = ^ , (Al5a) 

T>i = (A) • (A15b) 

In the case that O is strictly real on every configuration OD n is imaginary for odd n. In order to keep (O) real 
(OD n ) has to vanish for odd n and the preceding expansion simplifies to 

(0>„ = (O) (l + (-X> 2 + O2) 1? + {V\ - V A - V 2 2 + Oi) 

+ {-V\ + 22? 2 2? 4 -V 6 + V\0 2 - V 4 2 - V 2 0± + O e ) fi 6 ) + 0(p s ) , (A16) 

i.e. this formula is applicable to the correlation function in eq. 3. 

Because free energies are calculated from logarithms of correlation functions we give here the expansion coefficients 
of the logarithm of an observable O which can be obtained by inserting the above expansion into the expansion of 
the logarithm. For a generic, not necessarily real observable the expansion is 

In (O)^ = In (O) + Oin + (-V 2 - \o\ + 2 ) + Qo? - O1O2 + O3) ^ + Q©2 - ^4- 

x o\ + o\o 2 - \o\ - Oi0 3 + o 4 ) m 4 + (\o\ - 0\0 2 + OxOl + o\o 3 - o 2 o 3 



4 1 1 2 



3 



-O x O A + O b ) /x 5 + --2? 2 3 + 2? 2 2? 4 - As - -Of + 0\0 2 - -0\0\ + -0| - 0\0. 



+20i0 2 3 - ^3 + 01^4 - O2O4 - Oi0 5 + O e ) + ^(M 7 ) • (A17) 



For real observablcs this reduces to 



In (O)^ = In (O) + (-V 2 + 2 ) ^ 2 + Q©i - ©4 - ^0 2 2 + O4) M 4 



4 I + 2? 2 2? 4 - X» 6 + ^2 - 2 4 + O e I M B + O(m') • (A18) 
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